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Abstract.—Salamandra infraimmaculata is one of the prominent members of the Turkish herpetofauna, which 
is currently classified as Near Threatened by the IUCN. Although many studies of its morphology, ecology, 
and phylogeny have been conducted in recent years, many issues concerning its taxonomic structure and 
morphometry remain unresolved. In the present study, morphometric characters and color-pattern features 
of the specimens of Salamandra infraimmaculata captured from ivriz, Halkapinar, Konya Province are given 
in detail and compared with the data available in previous studies. In addition, the climatic niche preferences 
of this species were determined using ecological niche modelling, and three WorldClim bioclimatic variables 
were found to restrict the species presence: Minimum Temperature of Coldest Month (Bio6), Mean Temperature 
of Wettest Quarter (Bio8), and Precipitation of Warmest Quarter (Bio18). Even though this study revealed the 
potential habitat preferences of this species, clearly more detailed studies are needed to resolve the problematic 
taxonomical issues. 
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Introduction 


The genus Salamandra currently includes six recognized 
species: S. algira Bedriaga, 1883; 5. atra Laurenti, 1768; 
S. corsica Savi, 1838; S. infraimmaculata Martens, 1885; 
S. lanzai Nascetti, Andreone, Capula and Bullini, 1988; 
and S. salamandra Linnaeus, 1758 (Rodríguez et al. 
2017). One of them, the Near Eastern Fire Salamander 5. 
infraimmaculata, s a member ofthe Turkish herpetofauna 
(Baran etal. 2012). Salamandra infraimmaculata was first 
described as Salamandra maculosa var. infraimmaculata 
from Bcharré, Lebanon (Martens 1885). This species 
is now known from Turkey, Syria, Lebanon, northern 
Israel, northern Iraq, and western Iran and includes three 
subspecies: 5. i. infraimmaculata Martens, 1885, 5. i. 
orientalis Wolterstorff, 1925, and S. i. semenovi Nesterov, 
1917 (Joger and Steinfartz 1995; Steinfartz et al. 2000). 
In Turkey, these three subspecies are found in Hatay 
province (S. i. infraimmaculatay, in Adana, Mersin, and 
Malatya provinces of southern Turkey (S. i. orientalis), 
and in Bitlis and Erzincan provinces of eastern Turkey 
(S. i. semenovi) (Joger and Steinfartz 1995; Veith et al. 
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1998; Steinfartz et al. 2000). 

Specimens of this taxon from Turkey were previously 
treated as S. salamandra (Eiselt 1966; Schmidtler 
and Schmidtler 1970; Oz 1987; Oz and Arikan 1990; 
Arikan et al. 1990; Baran and Oz 1994). However, the 
Kemaliye (Erzincan) population was assigned to 5. s. 
infraimmaculata (Fachbach 1971), and the subspecies 
was later elevated to the full species level based on 
serum protein patterns (Gasser 1978). According to 
comparisons of blood-serum proteins among the various 
Salamandra populations by polyacrylamide disc gel 
electrophoresis, the populations in Turkey were classified 
as S. infraimmaculata (Joger and Steinfartz 1995). The 
study by Bóhme et al. (2013) on various attributes of 
S. infraimmaculata reported that S. i. semenovi has a 
typical scrolled pattern of yellow rings, semicircles, and 
similar pattern elements, while the color pattern of both 
nominates of form S. i. orientalis has larger broad solid 
flecks and small, yellow spots over the whole body except 
for the belly. In S. i. infraimmaculata, the yellow dots are 
large and extend over the whole body, except the belly. 
There are usually four yellow spots on the head, one on 
each parotoid and one above each eye (Rastegar-Pouyani 
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and Fizi 2006). Previous molecular phylogenetic studies 
identified six major monophyletic groups belonging to 
this genus, which were separated from each other by 5 
to 13 million years (Veith et al. 1998; Steinfartz et al. 
2000; Weisrock et al. 2006; Vences et al. 2014; Pyron and 
Wiens 2011; Rodríguez et al. 2017). 

One molecular study claimed that the populations 
occupying the Zagros Mountains (Iran) are genetically 
close to the southeastern populations of Turkey, and it has 
been suggested that some populations in southern Turkey 
can be considered as a new subspecies because the type 
locality of S. i. semenovi is in the other clade (Ahsani 
et al. 2019). Based on additional herpetological studies 
in Turkey, the known geographical distribution of S. 
infraimmaculata has recently been expanded (Coskun et 
al. 2013; Olgun et al. 2015; Sarikaya et al. 2017; Akman 
et al. 2018; Sami and Yildiz 2018; Yildiz et al. 2019), 
although the taxonomic state of the species division into 
three subspecies is still not fully clarified (Steinfartz et al. 
2000; Olgun et al. 2015; Ahsani et al. 2019). 

In this study, the morphological characters of S. 
infraimmaculata salamanders collected from a new 
locality in Turkey are described, which may shed some 
light on the taxonomy of this species. In addition, the 
climatic conditions affecting the distribution of 5. 
infraimmaculata were determined using ecological niche 
modelling. 


Materials and Methods 


Sampling. On 17 November 2019, one adult male (Fig. 
1А) and one adult female (Fig. 1B) of S. infraimmaculata 
were collected from Ivriz, Halkapinar, Konya Province 
(Fig. 2). The salamanders used in this study were 
found under stones during the daytime (between 1100 
and 1200 h) and collected by hand after a heavy rain. 
The localities where the specimens were collected are 
situated near small creeks, and their geographic positions 
were recorded with a GPS receiver (Garmin eTrex 30). 
The sex of the captured individuals was determined 
through secondary sexual characters, 1.е., males have a 
prominent (swollen) cloaca (Başoğlu and Özeti 1973). 
Morphometric measurements and the color and pattern 
characteristics of the specimens were recorded in the field 
and the specimens were released at the capture locations. 


Morphological characteristics. All morphological 
measurements were recorded using a digital caliper 
(Mitutoyo, Kawasaki, Japan) with an accuracy of 0.01 
mm. Measurements of body parts and their ratios follow 
previously published papers on salamanders (Oz 1987; 
Olgun et al. 2015), and are as follows: total body length 
(TBL), tip of snout to tip of tail; body length (LCP), 
length from snout to anterior end of cloaca opening; 
snout-vent length (SVL), tip of snout to posterior end of 
cloaca opening; tail length (TL), length from posterior end 
of cloaca opening to tip of tail; forelimb length (FLL); 
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hindlimb length (HLL); distance between fore- and hind 
limbs length (DFHL); head length (HL), distance from 
snout to gular fold; head width (HW); parotoid length (PL); 
parotoid width (PW); and distance between anterior of each 
parotoid (DAP). In addition, the ratios of DFHL/SVL, 
PL/HL, PW/HW, PW/PL, HW/HL, TL/SVL, HL/LCP, 
FLE/SVL. TL/SVL*100, HL/LCP*100, HW/HL*100, 
FLL/SVL*100, HLL/SVL*100, and TL/TBL*100 were 
calculated. Statistical analyses were carried out using 
STATISTICA 6.0 (StatSoft, Inc., Tulsa, Oklahoma, USA) 
to determine descriptive statistics for the measurements of 
the salamanders. Morphometric measurements were then 
compared to previously published data (Öz 1987; Öz and 
Arıkan 1992; Olgun et al. 2015). 


Ecological niche modelling. A total of 56 records 
were collected from the published literature and on- 
going fieldwork over the course of a few years (Fig. 2; 
Supplementary Table 1). The 56 records were thinned to 
44 localities using ‘spThin’ (Aiello-Lammens et al. 2015) 
in К (К Core Team, 2019), and the thinning distance 
was selected as 10 km. WorldClim bioclimatic variables 
(Bio1-19, Supplementary Table 2) were downloaded at 
0.5 arcmin resolution using the raster package (Hijmans 
2017) as environmental predictors. The minimum convex 
polygon was used as the background extent of the study 
region, and buffered by 0.5 degrees. Predictor rasters 
by the background extent and the random background 
points sampled were masked with values of 10,000. These 
processes were handled with the sp (Pebesma and Bivand 
2019) and rgeos (Bivand and Rundel 2019) packages in 
R. The occurrence records were partitioned by the block 
method (k = 4) using the ENMeval (Muscarella et al. 2014) 
package in R, and MaxEnt was run successfully with 
output evaluation results for 45 clamped models. Wallace, 
which is a flexible platform for reproducible modeling of 
species niches and distributions, was used for the complex 
workflows above (Kass et al. 2018). Later, according to the 
lowest AICc value (Supplementary Table 3), MaxEnt was 
run again with the following parameters: Linear, Hinge 
and Quadratic as features; regularization multiplier: 4; and 
number of replicates: 30. The Area Under the Receiver 
Operating Characteristic (ROC) Curves (AUC) value, 
averaged over the 30 replicated runs, was considered as 
an additional measure of model performance. Models 
with AUC = 0.5 indicate a performance equivalent to 
“random,” AUC > 0.7 indicates “useful” models, and 
AUC > 0.9 indicates models with “excellent” performance 
(Manel et al. 2001). 


Results 


Descriptive statistics of the metric measurements of 
the specimens are given in Table 1. The total length of 
specimens (TBL) was measured as 226.83 mm for the 
male and 155.85 mm for the female. The SVL of the male 
was 134.92 mm, while it was 99.25 mm for the female. 
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Distribution range expansion of Salamandra infraimmaculata 
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(A) male and (B) female. 


In the two specimens, the ground color of the upper 
side of the head is black with two larger solid yellow 
flecks (Fig. 1). In the male, these flecks are separated from 
each other, but in the female, they are in contact. While 
the female has small yellow spots on the side of the head 
on a black background, the male specimen 15 spotless. 
There is no pigmentation forming a faint dark dot on 
the parotoid glands. The ground color of the dorsum is 
blackish with a scattered pattern of large varying forms 
of yellow maculation that are sometimes in contact with 
each other. There are yellow spots on the side of the 
body, especially in the female. Although the gular region 
has maculation, the female has two obvious yellow flecks 
rather than dot-shaped staining, which is more common 
in males. In males, the ventral surface shows sparse 
dots consisting of small yellow spots, while females 
have purely black undersides. Both specimens show 
medium-sized yellow flecks on the black background of 
the extremities (Fig. 1). Comparative morphological data 
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on the populations of this species in different regions of 
Turkey are given in Table 2. 

The average test AUC for the replicate runs was 
0.863, and the standard deviation was 0.145. Among the 
bioclimatic variables, B106 showed the highest percentage 
contribution (51.9%), whereas B1018 and Bio8 had lower 
percentages of 20.7 and 12.1, respectively, and all other 
variables were under 10%. The optimal habitat for this 
species within the minimum convex polygon appeared 
to be across the middle Taurus and Amanos mountains in 
southern Anatolia (Fig. 3). 


Discussion 


Salamandra infraimmaculata is a polytypic species, 
which is currently recognized as three discrete 
taxonomic units: 5. i. infraimmaculata Martens, 1885 
(Turkish Hatay, Syria, Lebanon, and northern Israel), 5. 
i. orientalis Wolterstorff, 1925 (south and southeastern 


March 2022 | Volume 16 | Number 1 | e307 


Candan 


47 S.l semenovi 
S. i. orientalis 
CC S&Linfrimmaculata 


"4, 


e 


200 
= Kilometers 


Fig. 2. Distribution patterns 
with the new locality record. 
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Fig. 3. Predicted distribution of Salamandra infraimmaculata under current climatic conditions. Warm 
colors (red and yellow) show suitable habitats, whereas the blue color represents unsuitable habitats for 


S. infraimmaculata. 


Turkey), and S. i. semenovi Nesterov, 1917 (easternmost 
Turkey, western Iran, and northern Iraq) (Joger and 
Steinfartz 1995; Steinfartz et al. 2000; Bóhme et al. 
2013). However, the taxonomic status of its subspecies 
and their distributions are still unclear (Steinfartz et al. 
2000; Bóhme et al. 2013). 

The new locality record presented in this study 
extends the known range of S. infraimmaculata by about 
70 km to the northwest, as measured from the nearest 
previous locality of Findikipinari, Mersin. The new 
locality (Ivriz, Halkapinar, Konya) indicates that the 
distribution area of the species may extend throughout 
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the Eastern Taurus Mountains to the western direction. 
Several studies have reported data on the body sizes 
ofthe three subspecies. Schorn and Kwet (2010) reported 
that S. i. orientalis 1s smaller in size than the other two 
subspecies (S. i. infraimmaculata and S. i. semenovi). 
The maximum value for Israeli Fire Salamanders (5. 
i. infraimmaculata) was given as 316 mm by Eiselt 
(1958), while it was reported as 202 mm for Iraqi 
specimens corresponding to S. i. semenovi (Bóhme et 
al. 2013). Finally, within the distribution area of the 
S. i. orientalis, Bóhme et al. (2013) reported the total 
length of specimens from Findikpinari, Mersin, and Ilica 
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and blue shading respectively show the mean responses of the 30 replicate MaxEnt runs and the mean plus/minus one 


standard deviation. 


in Turkey as 200 mm and those from Kahramanmaras, 
Turkey as 255 mm. In the present study, the total body 
length was 226.83 mm for the male specimens. As can 
be seen from the values given above, the suggestion that 
adults of S. i. orientalis are smaller than the other two 
forms, as Bóhme et al. (2013) stated, is not correct. 
According to Bóhme et al. (2013), S. i. semenovi has 
the typical, scrolled pattern of yellow rings, semicircles, 
and similar pattern elements, while the color pattern of 
both S. i. orientalis and the nominotypical form consists 
of larger solid flecks. The color-pattern characteristics 
alone are not enough to distinguish the Fire Salamanders 
for taxonomically assigning individuals (Bóhme et al. 
2013). However, phylogenetic and phylogeographic 
information obtained from different molecular 
techniques can be combined with the knowledge of 
morphology and distribution, producing a more accurate 
taxonomic placement for the studied specimens. 
Considering the niche modelling of S. infraimmaculata 
in Syria, one study reported a positive relationship 
between fitness and the Precipitation of Coldest 
Quarter (Bogaerts et al. 2013). In this study, among 
the climatic variables used in the ecological niche 
modelling, Minimum Temperature of the Coldest 
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Month (Bio6), Mean Temperature of Wettest Quarter 
(B108), and Precipitation of Warmest Quarter (В1о18) 
are the most important. Generally, in southern 
Anatolia, S. infraimmaculata can survive in the low 
temperatures in winter months. In other words, the 
known occurrences for 5. infraimmaculata have 
minimum temperatures between -15 and 0 °C (Fig. 
4A). Similarly, S. infraimmaculata can tolerate mean 
temperatures between -8 and 0 °С in summer months 
(Fig. 4B). Additionally, the species 15 not found in 
localities with high precipitation, whereas it can adapt 
to low precipitation in the summer months (Fig. 4C). In 
particular, the temperature niches of S. infraimmaculata 
can explain local population abundances in the south of 
Anatolia (Bowler et al. 2015). 

There are still many uncertainties regarding the 
taxonomic structure of S. infraimmaculata, which 
makes it difficult to determine the distribution limits 
of the taxonomic groups within this species. In this 
study, a new locality record for the species is reported 
by revealing its potential habitat preferences. However, 
more detailed studies are needed to fully resolve 
the taxonomic uncertainties and to more thoroughly 
document the distribution preferences of the species. 
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Supplementary Table 1. Presence data of the Salamandra infraimmaculata localities used in the model. 


Locality 


Ivriz, Halkapinar, Konya 
Maden Village, Ulukisla, Nigde 
Camliyayla, Mersin 
Findikpinar, Mersin 

Gózne, Mersin 

Ömerli Village, Pozantı, Adana 
Omerli Village, Pozanti, Adana 
Kamışlı, Pozantı, Adana 
Akçatekir, Pozantı, Adana 
Belemedik, Pozantı, Adana 
Saimbeyli, Adana 

Feke, Adana 

Kozan, Adana 

Boztahta, Kozan, Adana 
Çulluuşağı, Kozan, Adana 
Camdere, Kozan, Adana 

Baraj, Kozan, Adana 

Yamanlı, Tufanbeyli, Adana 
Çamlıbel, Aladağ, Adana 
Belenköy, Adana 

Darıpınarı, Adana 

Bahçe, Osmaniye 

Kozaklı Village, İskenderun, Hatay 
İskenderun, Hatay 

Bekbele Plateau, İskenderun, Hatay 
Demen Plateau, İskenderun, Hatay 
Çardak Plateau, Hassa, Hatay 
Eğribucak, Hassa, Hatay 

Kapılı yolu, Dörtyol, Hatay 
Harbiye, Hatay 

Kozkalesi, Altınözü, Hatay 
Yeniköy, Samandağ, Hatay 
Arkıtça, Kırıkhan, Hatay 


Altınboğa Village, Andırın, Kahramanmaraş 


Ilıca Village, Kahramanmaraş 
Sugözü, Besni, Adıyaman 

Yağmurlu village, Gerger, Adıyaman 
Tut, Adıyaman 

Kaşlıca, Tut, Adıyaman 

Meydanköy, Gölbaşı, Adıyaman 
Harmanlı, Gölbaşı, Adıyaman 
Salihli Village, Kemaliye, Erzincan 
Kemaliye, Erzincan 


Yuva Village, Kemaliye Erzincan 
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Latitude 


37.430556 
37.448889 
37.183016 
36.917778 
37.005833 
37.528056 
37.526944 
3/52 
37.329722 
37.365833 
38.000833 
37.810833 
37.444722 
37.635556 
37.680278 
371653333 
374912222 
38.182307 
37.474429 
37.865833 
37.141883 
3142177222 
36.479444 
36.585556 
36.612550 
36.635556 
36.838611 
36.750833 
36.864444 
36.124722 
36.099980 
36.222778 
36.509167 
37.626096 
37.839878 
37.669669 
37.968325 
37.808611 
37.803813 
37.810277 
37.843888 
39.330269 
39.266111 
39.247924 


145 


Longitude 


34.181111 
34.630278 
34.591682 
34372223 
34.559722 
34.8575 
34.851111 
34.955556 
34.786389 
34.921944 
36.090556 
35.919167 
35.807778 
36.011388 
35.864167 
35.816389 
35.863889 
36.242130 
35.068174 
35,853277 
34.725625 
36.573055 
36.122778 
36.218333 
36.229809 
36.244167 
36.443611 
36.454166 
36.399167 
36.131944 
36.189077 
35.986389 
36.294999 
36.340641 
36.870419 
37.786397 
39.029137 
SE 
37-9179159 
37.644999 
37.758056 
38.502082 
38.4975 
38.509211 
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Reference 


This study 
Baran and Öz (1994) 
Baran and Öz (1994) 
Baran and Öz (1994) 

Olgun et al. (2015) 
Olgun et al. (2015) 
Olgun et al. (2015) 
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Baran and Öz (1994) 
Baran and Öz (1994) 
Sarıkaya et al. (2017) 
Sarıkaya et al. (2017) 
Sarıkaya et al. (2017) 
Sarıkaya et al. (2017) 
Sarıkaya et al. (2017) 
Sarıkaya et al. (2017) 
Sarıkaya et al. (2017) 
Böhme et al (2013) 
Özeti and Yılmaz (1994) 
Baran and Öz (1994) 
Olgun et al. (2015) 
Yıldız et al. (2019) 
Yıldız et al. (2019) 
Olgun et al. (2015) 
Yıldız et al. (2019) 
Yıldız et al. (2019) 
Yıldız et al. (2019) 
Yıldız et al. (2019) 
Yıldız et al. (2019) 
Yıldız et al. (2019) 
Olgun et al. (2015) 
Olgun et al. (2015) 
Sami and Yıldız (2018) 
Sami and Yıldız (2018) 
Sami and Yıldız (2018) 
Sami and Yıldız (2018) 
Sami and Yıldız (2018) 
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Supplementary Table 1 Continued. Presence data of the Salamandra infraimmaculata localities used in the model. 


Locality Latitude Longitude Reference 
Aslantepe, Malatya 38.381928 38.361328 Baran and Oz (1994) 
Gündüzbey Village, Malatya 38.278889 38.269444 Olgun et al. (2015) 
Eskihalfeti, Şanlıurfa 345989722 37.871666 Olgun et al. (2015) 
Tagar Stream, Cemisgezek, Tunceli 39.064444 38.905278 Olgun et al. (2015) 
Sütlüce, Tunceli 39.036944 39.165833 Olgun et al. (2015) 
Gelinodaları, Pülümür, Tunceli 39.470112 39.906898 Olgun et al. (2015) 
Soğuksu Village, Ergani, Diyarbakir 38.399089 39.656633 Olgun et al. (2015) 
Sağırkaya, Hizan, Bitlis 37.991108 42.569625 Akman et al. (2018) 
Alatepe Village, Ilıca, Bingöl 39.059281 40.756134 Çiçek et al. (2017) 
Derinçay Village, Karlıova Bingöl 39.132289 40.824913 Çiçek et al. (2017) 
Deliktaş, Bitlis 38.347686 42.041075 Akman et al. (2018) 
Kösepınarı, Osmaniye 37.595694 36.155118 Unpub. data (Serkan Gül) 


Supplementary Table 2. Definitions of the 19 WorldClim bioclimatic variables. 


Variable Definition 

BIO! Annual Mean Temperature 

BIO2 Mean Diurnal Range (Mean of monthly (max temp - min temp)) 
BIO3 Isothermality (BIO2/BIO7) (x100) 

BIO4 Temperature Seasonality (standard deviation x100) 
BIOS Max Temperature of Warmest Month 

BIO6 Min Temperature of Coldest Month 

BIO7 Temperature Annual Range (BIO5-BIO6) 
ВІО8 Mean Temperature of Wettest Quarter 

BIO9 Mean Temperature of Driest Quarter 
BIO10 Mean Temperature of Warmest Ouarter 
ВТО11 Mean Temperature of Coldest Quarter 
BIO12 Annual Precipitation 

BIO13 Precipitation of Wettest Month 

BIO14 Precipitation of Driest Month 

BIOI5 Precipitation Seasonality (Coefficient of Variation) 
BIO16 Precipitation of Wettest Quarter 

BIO17 Precipitation of Driest Quarter 

BIO18 Precipitation of Warmest Quarter 

BIO19 Precipitation of Coldest Quarter 
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Supplementary Table 3. Performance of the 48 models created during the evaluation process for Salamandra infraimmaculata. The 
gray-shaded entry shows the settings chosen for the model. FC: Feature Classes, RM: Regularization Multiplier, AICc: corrected 
Akaike Information Criterion. 


FC RM AICc FC RM AICc FC RM AICc 
L 0.5 1091.86 Г 2 1094.35 Г 3.5 1095.11 
LQ 0.5 1076.99 LQ 2 1053.38 LQ 3.5 1055.45 
H 0.5 NA H 2 1078.45 H 3.5 1065.6 
LOH 0.5 NA LOH 2 1045.37 LQH 3.5 1045.81 
LOHP 0.5 NA LOHP 2 1042.53 LOHP 3.5 1051.03 
L 1 1087.92 L 2.5 1090.27 L 4 1096.62 
LO 1 1059.28 LO 25 1053.95 LQ 4 1056.35 
H 1 1383.18 H 2.5 1067.06 H 4 1055.19 
LOH 1 1289.72 LOH 2.5 1049.89 ^. LOH 4 1040.85 
LOHP 1 1332.03 LOHP Do 1049.21 LOHP 4 1065.72 
L 1.5 1088.85 L 3 1091.14 L 4.5 1101.69 
LQ 1.5 1060.73 LQ 3 1054.65 LQ 4.5 1057.37 
H 1.5 1087.01 H 3 1061.92 H 4.5 1058.83 
LOH 1.5 1057.83 LOH 3 1051.48 LOH 4.5 1042.09 
LOHP 1.5 1069.73 LOHP 3 1048.21 LOHP 4.5 1051.51 
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